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Abstract 

The present study is aimed at investigating the validity of the generalized second law (GSL) of thermodynamics in 
f(H) gravity. Choosing f(R) = 11 + §72/* + ^1Z~ V with §, /x, f , v > (following the study of Nojiri and Odintsov 
in 2003), we have computed the time derivatives of total entropy for various choices of scale factor pertaining to 
emergent, intermediate, and logamediate scenarios of the universe. We have taken into account the radii of Hubble, 
apparent, particle, and event horizons while computing the time derivatives of entropy under various situations being 
considered. After analyzing through the plots of time derivative of total entropy against cosmic time, it is observed 
that the derivative always stays at positive level. This indicates the validity of GSL of thermodynamics in the f(7L) 
gravity irrespective of the choices of scale factor and enveloping horizon. 
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Background 

Accelerated expansion of the universe is well documented 
in literature (detailed discussion is available in [1] and ref- 
erences therein). Approximately 76% of the energy content 
of the universe is not dark or luminous matter, but it is 
instead a mysterious form of dark energy that is exotic, 
invisible, and unclustered [2] . In order to explain the accel- 
erated expansion of the universe, three main classes of 
model exist [3]: (1) a cosmological constant A, (2) dark 
energy, and (3) modified gravity. The last class of mod- 
els known as extended theories of gravity corresponds to 
the modification of the action of the gravitational fields 
[3]. These theories are based on the idea of an exten- 
sion of the Einstein-Hilbert action by adding higher order 
curvature invariants. Modified gravity theories have been 
reviewed in [4-6]. Nojiri and Odintsov [7] suggested f(7Z) 
gravity characterized by the presence of effective cosmo- 
logical constant epochs in such a way that early time 
inflation and late time cosmic acceleration are mutually 
unified within a single model. In another work, Nojiri 
and Odintsov [8] proposed another class of modified 
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f(TZ) to unify 7Z m inflation with A -cold dark matter era. 
Chattopadhyay and Debnath [9] considered f(TZ) gravity 
in an universe characterized by a special form scale factor 
known as emergent scenario and observe and concluded 
that the EoS parameter behaves like quintessence in this 
situation. 

The relevance of studying the thermodynamics of the 
universe was precisely mentioned in [10] that investigated 
the validity of generalized second law (GSL) in f(1Z) grav- 
ity. In the present work are going to investigate the GSL 
of thermodynamics for various choices of the scale factor 
a. The choices are named as emergent', 'intermediate' and 
logamediate! The physical aspects behind such choices 
are well documented in the literature [11-13]. However, 
for the sake of convenience we shall give a brief overview 
in a subsequent section. In this work the thermodynamic 
consequences of the universe for the said choices of scale 
factor would be examined in a modified gravity theory 
named as/(7£) gravity that has gained immense interest 
in recent times. In the present work we have extended 
the study of [9] to the investigation of the generalized 
second law (GSL) of thermodynamics in emergent sce- 
nario with the universes enveloped by Hubble, apparent, 
particle, and event horizons, respectively. In the remain- 
ing part of the paper, the radii of the said horizons are 
denoted by R\\, Ra> Rv> and Re, respectively. The validity 
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of the GSL implies that the sum of the time derivatives of 
the internal entropy and entropy on the horizon is non- 
negative. Hence, the primary objective of this work is to 
discern whether S total = ^internal + ^horizon > holds 
for the situations under consideration. The GSL would be 
investigated based on the first law of thermodynamics. 
Relevance of the laws of thermodynamics in cosmology 
was discussed by [14,15]. The validity of GSL in various 
DE candidates and their interactions has been discussed 
in several papers [16-22]. The works on the validity of the 
GSL in modified gravity theories include [23-27]. Bamba 
and Geng in [28] studied thermodynamics of the apparent 
horizon in f(7l) gravity, and it was demonstrated that an 
f(JZ) gravity can realize a crossing of the phantom divide 
and can satisfy the second law of thermodynamics in the 
effective phantom phase as well as non-phantom one. 
Another work [25] studied the thermodynamic behav- 
ior of field equations for f(lZ) gravity. Sadjadi and Jamil 
[29] applied logarithmic correction to the usual form of 
entropy and investigated the conditions that the presence 
of such modified terms in the entropy puts on other physi- 
cal parameters the system such as the temperature of dark 
energy via requiring the validity of GSL. In another recent 
work, Sadjadi [30] considered spatially flat Friedmann- 
Robertson-Walker (FRW) universe in the framework of 
the modified Gauss-Bonnet gravity and obtained the con- 
ditions required for validity of generalized second law 
(GSL). 

In the present work, we have taken different choices 
for the scale factor and examined whether the GSL holds 
for those choices. Details are presented in the subsequent 
sections. In a very recent work, Debnath et al. [31] investi- 
gated the validity of GSL for various choices of scale factor 
in fractional action cosmology by using, as well as without 
using, the first law of thermodynamics and considering 
different enveloping horizons. The present work is devi- 
ated from the study of [31] in the sense that instead of 
considering fractional action cosmology, we have chosen 
f(1Z), a modified gravity, to investigate the validity of GSL. 

Methodology 

The generalized second law 

In this section we are going to examine whether the gen- 
eralized second law (GSL) will hold for various choices of 
scale factor and on various horizons under f(JZ) gravity. 
The basic necessity for the validity of GSL is that the time 
derivative of the total entropy Stotal = 5h + S > 0, where 
S indicates the time derivative of normal entropy and 5h 
indicates the horizon entropy [22]. 

We consider the general class of FRW models with 
scalefactor a(t) [32], 



ds 2 = dt 2 + a 2 (t) [ X 2 + Sl( X )(d0 2 + sin 2 <9^ 2 )] (1) 



where S#(x) = sm X> X> sinhx are for closed, flat, and 
open models, respectively; and a(t) is the scale factor. 
For non-flat FRW universe, the field equations are 



H z + — = 



SnG 



H — 



-4nG(p+p) 



(2) 



(3) 



The first law of thermodynamics (Clausius relation) [33] 
on the horizon is denned as TxdSx = 8Q = —dE\, where 
dE is the amount of the energy flow through the cosmo- 
logical horizon. From the unified first law, reference [34] 
derived that during the time internal dt, one can obtain 
the amount of energy crossing the cosmological horizon 



-dE x = T x dS x = A c R x H(p + p)dt, 



(4) 



where Tx and Rx are the temperature and radius of the 
horizons under consideration in the equilibrium ther- 
modynamics. The area A c of the cosmological horizon 
generalized to curved space is [32] 

A c = 47ta 2 (t)Sl(Xc) (5) 

Subsequently, from Equations 4 and 5 the time deriva- 
tive of the entropy on the horizon can be derived as 



Sx 



47TRxa 2 (t)S 2 k (xc)H 
Tx 



(P+P) 



Using Equation 3 in Equation 6 we get 

Rxa 2 (t)S 2 k ( X c)H 



Sx = 



GTx 



(*-*) 



(6) 



(7) 



Time derivative of entropy inside the horizon is given 
by [35] 



TSix = (p + p)dV c + V c dp, 



(8) 



where volume V c is within the cosmological horizon. 
Using the field equation we get 



Six 
where [32] 

Vc 



-C-7 2 ) 



(Vc - 3HV C ) 



27ra 3 (t)(xc — sin Xc cos Xc) closed, 



4^3 



a\t)xc flat, 
2na 3 (t)(-Xc + sinh Xc cosh Xc) open. 



(9) 



(10) 



In the present work, we are going to consider four 
cosmological horizons, namely, (a) apparent, (b) Hubble, 
(c) event, and (d) particle horizons. For the above hori- 
zons, the radii are expressed as R\, Ru, R^ and i?p, 
respectively. The above horizons are given by [35] 

£a = , 1 (11) 



+ ^ 
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Here, H = | is the Hubble parameter. For k = (i.e., flat 
universe) we get the radius of the Hubble horizon. Thus, 
the radius of the Hubble horizon is Rh = jj. 



Re 



Rj> 



f°° dt 
= a I — 
St a 

f l dt 
= a I — 
Jo a 



R A = —HR K \H - — J ; R H = - 
R E = HR E - 1 ; R P = HRp + 1 ; 



H 

'Tf 2 ' 



(12) 



(13) 



(14) 



Our target is to investigate whether Sx + Six > holds 
for the cases to be considered. 

Basic equations of f (11) gravity 

The action oif(JZ) gravity is given by [28] 



j d 4 x^ 



2k 2 



(15) 



where g is the determinant of the metric tensor g^ vi 
^matter is the matter Lagrangian and k 2 = 8ttG. The 
f(7Z) is a nonlinear function of the Ricci curvature 1Z that 
incorporates corrections to the Einstein-Hilbert action 
which is instead described by a linear function/(7£). The 
gravitational field equations in this theory are [28] 



H 2 + — = 



H 



SfOZ) 



(P + Pc) 



(p +p + p c +Pc) 



(16) 



(17) 



where p c and p c can be regarded as the energy density and 
pressure generated due to the difference oif(JZ) gravity 
from general relativity given by [28] 



Pc 



Snf _ 



f-Rf 



- mf'k 



Pc = 



f-nf 



+f"n +f"n 2 + 6f'n 



(18) 



(19) 



where the scalar tensor K = -6 (k + 2H 2 + ^j. 
The choices of scale factor a(t) 

In this paper we have considered three forms of the scale 
factor a(t) in the f(TZ) gravity to investigate the valid- 
ity of the GSL of thermodynamics. The three choices, in 



literature, dubbed as emergent, intermediate, and logame- 
diate respectively, are given by 

1. Emergent, a(t) = A (ri + e 8 *)" with 
A>0, B>0, r\ > 0, n>l [11]. 

2. Intermediate, ait) = expCB^ 3 ) with 
B > 0; < < 1 [12]. 

3. Logamediate, a(t) = exp(A(\nt) a ) with 
Aol > 0, a > 1 [13]. 

For the above choices of scale factor, the forms of the 
Hubble parameter H are the following 



H = 



Bne 31 
rj + eP*' 



H = B/3r 1 +P; H = 



Act (In t) 



-l+a 



(20) 



It is clear from the above equations that we are first 
choosing various forms of scale factor and subsequently 
investigating the GSL in the corresponding scenarios. This 
reverse' way of investigations had earlier been used exten- 
sively by Ellis and Madsen [36], who chose various forms 
of scale factor and then found out the other variables from 
the field equations. 

We choose the function f(lZ) as [37] 



f(K) = n + $w L + zn- 



(21) 



We obtain the Ricci scalar 1Z for the above three choices 
of scale factor leading to the forms of H obtained in (11). 
Subsequently we obtain f(lZ) for all of the above choices 
as functions of time t. The radii of the various enveloping 
horizons of the universe are given below. 

The radius of the apparent horizon is given by 



1 



H 2 + 



If we use k = 0, then we get the radius of the Hubble 
horizon R H = j^. The radii of the particle R]> and the event 
R E horizons are given by 



[ f dt , f 00 dt 

Rp = a I — and R E = a — . 
Jo a Jo a 



Discussions on the above radii of different horizons are 
available in [35]. Using the above forms of scale factors, 
the Ricci scalar 7Z is reconstructed as follows: 

For emergent scenario: 



71 = 6 



B 2 e 2Bt n 



+ 



2B 2 e 2Bt n 2 



+ 



(e Bt + rj) 2 (e Bt + T]) 2 

k(^ + r])- 2n 
A 2 



+ 



B 2 e Bt n 

e Bt _j_ ^ 



(22) 
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Figure 1 S t otai against cosmic time f for Hubble horizon in 
emergent scenario. 
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Figure 3 S tota i against cosmic time r for Hubble horizon in 
logamediate scenario. 



For intermediate scenario: 



U = 6 \e~ 2Bt ^k + Br 2 +P(-1 + p)p + 2£ 2 r 2+2/ ^ 2 ] 

(23) 



For logamediate scenario: 



K = 6 



-2A(\nt) c 



k + 



A(-l + a)a(\nt)- 2 + a 



Aa(\nt)~ l+a 2A 2 a 2 (\nt)- 2+2a 



(24) 



t 2 



t 2 



Now we have discussed the validity of the GSL of ther- 
modynamics with the various choices of scale factor by 
obtaining the time derivatives of total entropy from for 
the universe enveloped by the different horizons and then 
plotted the time derivatives of the total entropy against 
cosmic time t to get the following twelve graphs shown in 
the Figures 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 12 [Figures 1, 
2, and 3 (for Hubble horizon), Figures 4, 5, and 6 (for 
apparent horizon), Figure 7, 8, and 9 (for particle hori- 
zon), and Figure 10, 11, and 12 (for event horizon)]. In 



all the plots we find that ^total is staying in the positive 
level. This indicates the validity of GSL of thermodynam- 
ics in all scenarios of the universe which are enveloped 
by Hubble, apparent, particle, and event horizons. In all 
the figures the red, green, and blue lines correspond to 
k = — 1, +1, and k = 0, respectively. 

Results and discussion 

In the present work we have investigated the validity of 
generalized second law of thermodynamics in an universe 
enveloped by Hubble, apparent, particle and event hori- 
zons. Instead of considering FRW universe governed by 
Einstein gravity, we have considered a modified gravity in 
the form oifilZ) gravity. We have chosen the scale factors 
in three forms corresponding to emergent, intermediate, 
and logamediate scenarios. While investigating the valid- 
ity of GSL of thermodynamics, we have not taken into 
account the first law of thermodynamics. For the purpose 
of the investigation of the validity of GSL, we have com- 
puted the entropy on the horizon, as well as inside the 
horizon, in the all twelve cases under consideration. We 
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Figure 2 5 to tai against cosmic time f for Hubble horizon in the 
intermediate scenario. 
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Figure 4 5 tota | against cosmic time t for apparent horizon in the 
emergent scenario. 
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t 

Figure 5 S tota | against cosmic time f for apparent horizon in the 
intermediate scenario. 
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Figure 7 S tota | against cosmic time f for particle horizon in the 
emergent scenario. 



have kept the curvature of the universe under considera- 
tion. In all the possible three cases, we have examined the 
GSL for flat (k = 0), open (k = — 1), and closed (k = 1) 
universes. We have plotted the time derivative of the total 
entropy against the cosmic time t, in all of the cases 
under consideration. 

In Figures 1, 2, and 3 we have considered three choices 
of scale factors in a universe enveloped by Hubble horizon 
and characterized by f(7Z) gravity. In all of the three cases, 
^total remains at the positive level and exhibit decaying 
behavior with the passage of cosmic time t. This indicates 
validity of GSL of thermodynamics in a universe charac- 
terized by f(JZ) gravity and enveloped by Hubble horizon. 
Moreover, this holds irrespective of the curvature of the 
universe. It is further noted that for the intermediate and 
logamediate scenarios, the rate of decay of Stotal is faster 
than in the case of emergent scenarios. It is further noted 
that this decaying behavior is significantly influenced by 
the curvature of the universe, in the case of logamedi- 
ate scenario. From Figure 3 we find that, in logamediate 
scenario, ^total falls ver T sharply in case of flat universe 




t 

Figure 6 S tota | against cosmic time f for apparent horizon in the 
logamediate scenario. 



(/c = 0). However, this rate is much lesser in open (k = — 1) 
and closed (k = 1) universes. 

In Figures 4, 5, and 6 we have considered apparent 
horizon. Although «S t otal stays at the positive level in all 
these three cases, the nature of its decay with cosmic 
time t has varied with the choice of scale factor. It has 
been observed that, in the case of logamediate scenario 
(Figure 6), 5 to tal nas fallen very sharply irrespective of the 
curvature. Whereas, in the the case of emergent scenario 
(Figure 4), the rate of change of S tot2L \ is much lesser. In this 
case the decaying of S to tal is very slow for flat (k = 0). In 
the case intermediate scenario (Figure 5), the decaying of 
<S to tal is not significantly influenced by the curvature of the 
universe. 

Figures 7, 8, and 9 confirm the validity of GSL of ther- 
modynamics in a universe characterized by f(7Z) gravity 
and enveloped by particle horizon. Here .Stotal has not 
shown any significant dependence on the curvature of 
the universe. However, 5 to tal behavior exhibited significant 
changes in different scenarios of the universe. In the case 
of emergent scenario, it is increasing with cosmic time t; 
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Figure 8 5 to tai against cosmic time t for particle horizon in the 
intermediate scenario. 
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Figure 9 S tota | against cosmic time f for particle horizon in the 
logamediate scenario. 



but in the case of intermediate scenario, it is decaying with 
cosmic time t. Although in the case of logamediate sce- 
nario (Figure 9), «S to tal behaves differently from the other 
scenarios. In Figure 9 we can see that 5 to tal is decaying 
with cosmic time t after increasing up to a certain period 
of time. 

In Figures 10, 11, and 12 we have plotted the time 
derivatives of total entropy for the universe in the emer- 
gent, intermediate, and logamediate scenarios, respec- 
tively, for the universe enveloped by the event horizon. 
These figures reveal that in f(TZ) gravity the GSL of 
thermodynamics is valid for all the scenarios under con- 
sideration when we are assuming event horizon as the 
enveloping horizon of the universe. 

Therefore, the rigorous study reported above reveals 
the validity of the GSL of thermodynamics in a universe 
governed by f(7Z) gravity. Irrespective of the choice of 
scale factor, enveloping horizon and curvature of the uni- 
verse, the time derivative of total entropy stays at the 
positive level. 




t 

Figure 1 S t otai against cosmic time f for event horizon in the 
emergent scenario. 
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Figure 1 1 S t otai against cosmic time f for event horizon in the 


intermediate scenario. 



Conclusions 

While concluding the paper we now present a compar- 
ative analysis of the outcomes of the present paper with 
the existing works in this direction. In [28] the validity 
of GSL was investigated iovfiJZ) gravity on the apparent 
horizon, and it was shown that the GSL can be satisfied 
in both phantom and non-phantom phases of the uni- 
verse. The present study deviates from the said study in 
the respect for choosing a model of f(1Z) and considering 
various forms of the scale factor available in the litera- 
ture. Moreover, here we have not confined ourselves to 
the apparent horizon only. We have also considered the 
other enveloping horizons like Hubble, particle, and event 
horizons. In all of our cases under consideration, the GSL 
of thermodynamics has been found to be satisfied. This 
study has adopted an approach similar to that of Debnath 
et al. [31]. However, the present study deviates from the 
said work in its choice of field equations. In the present 
work f(1Z) has been considered, whereas in [31] the frac- 
tional action cosmology was considered. The conclusions 
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Figure 1 2 5 tota | against cosmic time f for event horizon in the 
logamediate scenario. 
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are deviated accordingly. In [35] logamediate and inter- 
mediate scenarios were considered in Einstein gravity and 
the validity of GSL was investigated. In the present work 
we observe that like in [35], the GSL is satisfied for both 
of the scenarios in spite of considering f(7Z) instead of 
Einstein gravity. 

In an old paper Barrow [38] looked at the validity of 
GSL in the context of viscous universes, where there is 
more diverse behavior and also the possibility for the cre- 
ation of a dissipative structure. A study, highly relevant 
to be mentioned here, was done by Sadjadi [10], where 
the GSL was investigated inf(7l) with the scale factor 
a(t) = ao(t s — t)~ n with n > and derived some con- 
ditions for the validity of GSL, supposing the temperature 
to be proportional to the Gibbons-Hawking temperature 
and taking future event horizon as the horizon of the uni- 
verse. In [10] f(K) = all m and/OK) = pK + all m were 
considered with a, /3, and m e ?fi and finally, some con- 
ditions on a, /3, and m were derived analytically for the 
validity of GSL. However, in the present paper we have 
not worked on deriving conditions on the model parame- 
ters. Rather, we have considered f(K) =R + + £R~ V 
with §, /x, £ , and v > following [37] and studied various 
forms of scale factor a(i) taken from the existing litera- 
ture. Instead of deriving conditions for model parameters, 
we have created plots for the time derivatives of the total 
entropy «S t otal and examined whether it is staying at posi- 
tive levels or not. As future study we propose to work for 
generalization of the results and to derive restrictions on 
the model parameters required for GSL to be validated in 
various scenarios. 
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